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We propose a novel kinetic exchange model differing from previous ones in two main aspects. 

First, the basic dynamics is modified in order to represent economies where immediate wealth 
exchanges are carried out, instead of reshufflings or uni-directional movements of wealth. Such 
dynamics produces wealth distributions that describe more faithfully real data at small values of 
wealth. Secondly, a general probabilistic trading criterion is introduced, so that two economic units 
can decide independently whether to trade or not depending on their profit. It is found that the 
type of the equilibrium wealth distribution is the same for a large class of trading criteria formulated 
in a symmetrical way with respect to the two interacting units. This establishes unexpected links 
between and provides a microscopic foundations of various kinetic exchange models in which the 
existence of a saving propensity is postulated. We also study the generalized heterogeneous version 
of the model in which units use different trading criteria and show that suitable sets of diversified 
parameter values with a moderate level of heterogeneity can reproduce realistic wealth distributions 
with a Pareto power law. 


I. INTRODUCTION 


Kinetic exchange models provide a minimal descrip¬ 
tion of wealth exchange between economic units (rep¬ 
resenting, e.g., individuals, families, or companies) in a 
way formally similar to that in which energy is trans¬ 
ferred between molecules of a fluid due to collisions 00 - 
Such models were introduced independently in different 
fields such as social sciences 0 - 01 , economics 0 - 0 , and 
physics (Ml. John Angle 0 |4] originally introduced 
this type of models basing them on the surplus theory, 
with the goal of describing the origin of wealth inequal¬ 
ities. Compared with other agent-based models of fi¬ 
nancial markets [ 10 , the structure of kinetic exchange 
models is extremely simple: they only describe wealth 
flows between economic units without considering (ex¬ 
plicitly) other elements of a market dynamics. However, 
even so they predict quite realistic shapes of wealth dis¬ 
tributions jig . 

In kinetic exchange models it is assumed that the total 
wealth of the system is conserved, following in turn from 
wealth conservation during each unit-unit interaction. 
The dynamics can be formulated through the following 
update rule at a generic time iteration t defining the 
wealth exchange between two units j and k chosen ran¬ 
domly among the N units of the system (j, k = 1,..., iV), 


Xj = Xj + Ax jk , 
x k = x k - A x jk ■ 

Here Xj and x k (Xj and x k ) are the wealths of the trading 
units j and k before (after) the exchange and A Xj k = 
—A x k j is the exchanged amount of wealth. Without loss 
of generality, in Eqs. m the plus (minus) sign has been 
chosen for j (fc); who wins or loses depends on the sign of 
A Xj k . The conservation of the total wealth of the system 
leads to an equilibrium wealth distribution that coincides 


formally with the Gibbs energy distribution. 

As discussed in Ref. 0 ] , most of kinetic exchange mod¬ 
els can be formulated in a unified way, where the amount 
of exchanged wealth in Eqs. m is expressed as 

Axj k — uj j k x k aj k j Xj . (2) 

Here ujj k and Cj k j are two suitable stochastic variables 
in (0, l) to be extracted at each trade, representing the 
fractions of wealth of unit k and j, respectively, involved 
in the trade. 

Various aspects of kinetic exchange models have been 
criticized. A first objection concerns wealth conservation. 
Namely, it is pointed out that the assumption of a perfect 
conservation of wealth is incompatible with the fact that 
in real systems wealth is not conserved, a most relevant 
reason being the existence of production-consumption 
processes. In fact, kinetic exchange models can be in¬ 
terpreted as models with homogeneous production and 
consumption. To show this, one can use a modified form 
of a kinetic exchange model where production and con¬ 
sumption have been explicitly added to the dynamics. It 
is convenient to start from Eqs. m and © and write the 
total wealth variation of a generic unit j during a time 
step Jf, defined as the time interval corresponding to one 
Monte Carlo swap over all the units k (k ^ j), 

x'a — Xj 

~ y [ JjkXk ~ JkjXj] + (p — c)Xj . (3) 

Here the homogeneous production and consumption 
terms have been added, with production and consump¬ 
tion rates p and c, respectively; Jj k = u>j k /5t (and 
analogously for Jkj)- In the continuous time limit, 
( x'j — Xj)/St ~ dxj(t)/dt , one can introduce the auxil¬ 
iary variables 

Xj(t) = Xj(t) exp J dt\p(t) - c(t)] , (4) 
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where a possible time dependence of p and c has been 
taken into account. Equations m written in the new 
variables Xj do not contain anymore the production and 
consumptions terms and turn into the equations of the 
corresponding model without production and consump¬ 
tion, in which wealth conservation holds. However, this 
is not valid in the case of strong heterogeneity, in which 
the consumption and production processes of each unit 
have to be taken into account separately. 

Another critic toward kinetic exchange models claims 
that the type of trade dynamics they model can hardly 
resemble an actual economic trade. In fact, traditional 
approaches of economics assume that the decisions to 
carry out a trade are taken by rational agents or by eco¬ 
nomic agents with bounded rationality on the base of the 
total or partial information available about the system. 
Instead, kinetic exchange models may resemble at first 
sight more a hazard game, due to their random dynam¬ 
ics mm. It is clear that due to their statistical nature 
kinetic exchange models certainly do not provide a direct 
picture of the economic trading or the exchange activity 
that they are supposed to describe. So far a satisfac¬ 
tory microscopic picture and justifications of the models 
are missing. Their justification is pragmatic in nature, 
in that it relies a posteriori on the successful prediction 
of realistic shapes of wealth distributions. A notable ex¬ 
ception is the investigation reported in Ref. 0 , where 
a direct link with micro-economics is suggested. 

The main goal of the present paper is to make a step 
toward a microscopic foundation of kinetic wealth ex¬ 
change models. This is carried out along a twofold path. 
As a first element of the proposed reformulation, a novel 
“immediate-exchange” dynamics is introduced in Sect. HD 
in order to represent an actual exchange of goods rather 
than a random reshuffling or a uni-directional flow of the 
wealths of the interacting units, characteristic of other ki¬ 
netic exchange models. This novel dynamics is shown to 
produce equilibrium wealth distributions that better de¬ 
scribes empirical data also in the very small wealth range. 
Secondly, in Sect. Mil we introduce a probabilistic crite¬ 
rion that is used by each economic unit in order to decide 
whether to carry out an exchange or not. The exchange 
is carried out only if both units accept to do it. Such 
an acceptance criterion describes the decision process of 
a single unit during each interaction and can be suitably 
customized to represent the type and amount of partial 
information available to a unit during the trade. In the 
present work we follow a microscopic approach consid¬ 
ering acceptance criteria only based on the information 
concerning the ongoing trade, that is directly available 
to the two interacting economic units. To the best of our 
knowledge a probabilistic approach based on the (par¬ 
tial) information available about the quality and price of 
a product has been introduced in micro-economic mod¬ 
els in Refs. 0,0, but the probabilistic criteria used in 
many-agents models usually follow an approach based on 
a utility function 0]. 

As an important check versus real data, the heteroge¬ 


neous version of the new model is studied in Sect. m 
It is shown that, analogously to e.g. the Chakrabarti- 
Chakraborti kinetic exchange model, there exist suit¬ 
able sets of diversified parameters reproducing realis¬ 
tic wealth distributions on all scales including the large 
wealth range where the Pareto power-law is observed. 

Results and possible future lines of research are dis¬ 
cussed in the Conclusion. 


II. IMMEDIATE-EXCHANGE MODEL 
A. Formulation of the immediate-exchange model 

In this section a kinetic exchange model that provides a 
description of a market where encounters between units 
are accompanied by immediate exchanges, is proposed. 
By “immediate exchange” it is meant for example the 
type of interchanges that characterize barter where goods 
are directly exchanged without using a medium of ex¬ 
change, or market economies where goods are exchanged 
with money immediately or according to some agreed 
time schedule. Immediate exchanges are to be contrasted 
with the “delayed exchanges” characterizing e.g. gift 
economies [0 , where for cultural reasons valuables are 
given without an explicit agreement for immediate or fu¬ 
ture rewards. In the latter case one can talk about unidi¬ 
rectional trades. Unidirectional trades can also be used 
to describe, e.g., insurance business. 

For the sake of clarity, we start by considering a 
barter model, in which one can think of each unit i 
(i = 1,2,..., N) as having some items that he is will¬ 
ing to change for something else at some point; the to¬ 
tal value of these items is x.j. The model is evolved in 
time by extracting randomly two units j and k at ev¬ 
ery time iteration. The two units interchange something 
with value ejXj and tkXk, respectively, where €j and e*, 
are two independent uniform random numbers in the in¬ 
terval (0,1], different at each iteration. The random na¬ 
ture of these quantities describes the situation in which 
each time when somebody is trying to do a trade it can 
be of a different object among all the things he owns and 
is willing to exchange for something else. 

Notice that the model does not distinguish between 
different types of goods but focuses only on the corre¬ 
sponding amounts of wealth exchanged between units; 
therefore, one of the goods can be possibly understood 
as the currency in use. Thus, the model can be rein¬ 
terpreted as the one of a currency-based market. In the 
case of a currency-based trade, where an object of value 
Xj is exchanged with some money Xk 1 the random num¬ 
bers describe intrinsic fluctuations of prices and currency 
values. 

The dynamics of the proposed model can be defined 
by the following equations: 

Xj = (1 - €j)xj + e k x k , 
x k = (l- e k )x k + CjXj , 


( 5 ) 
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which can be rewritten in the same form of Eqs. 0 with 

A Xj k = £kXk - tjXj . (6) 

Because ej , e k > 0 then from Eqs. 0 together with © 
it is clear that a situation where Xi = X (i.e., a unit i 
owns all the wealth of the system) or Xi = 0 (i.e., a unit 
i is totally poor) is not possible at any moment of time. 
In fact, numerical simulations show that the equilibrium 
distribution of wealth is a F -distribution, 


f*A x ) = Jrr(P x ) a 1 exp(-/3x), (7) 

r(a) 

with a = /? = 2. In the following, using the fact that 
(, x) = a//3 and the average wealth is constant and always 
set here as ( x) = 1 due to the initial conditions x, = 1, 
one has a = (3 and it is convenient to use the simplified 
one-parameter form f a (x) = [f a ,a/{x)( x )](x)= i, i.e., 


f a (x) 


OL 1 


orx 


r(a) 


exp(— ax). 


( 8 ) 


The equilibrium wealth distribution corresponding to the 
value of the shape parameter a = 2 found from the 
fitting for the immediate-exchange model is / 2 (x) = 
4xexp(—2x). Thus, the wealth distribution is zero at 
zero wealth, / 2 (0) = 0, and decays exponentially at large 
x. 

It is to be noticed that the value of a ~ 2 is in the 
range (2, 2.5) of the values characterizing real distribu¬ 
tions, such as the household incomes analyzed by Salem 
and Mount sm. This value has been obtained using 
the plain immediate-exchange dynamics without further 
hypotheses, as for instance the introduction of an explicit 
saving propensity constraining the amount of wealth en¬ 
tering an exchange, as in the Chakraborti-Chakrabarti 
model. 

As a technical remark, the considerations above are 
valid for random numbers e, uniformly distributed in 
(0,1). Employing random numbers with different dis¬ 
tributions will produce accordingly a modified shape of 
the equilibrium wealth distribution. Details about the 
numerical algorithms used in the simulations of kinetic 
exchange models can be found in Ref. [§| . 

We also notice that the dynamics described by Eqs. 0 
is in the very spirit of kinetic exchange models, reflecting 
even more closely than other kinetic exchange models the 
statistical similarity of wealth flows to inter-molecular en¬ 
ergy exchanges. In fact, according to kinetic theory, the 
energy exchanged in a collision between two molecules j 
and k has the form given by Eq. ©• The variables ej and 
e k depend then on the initial directions of the molecular 
velocities and are to be considered independent random 
numbers in the hypothesis of molecular chaos. For fur¬ 
ther details see Ref. [25| . 


B. Immediate-exchange versus unidirectional and 
reshuffling models 


The immediate-exchange model introduced above 
is formally very similar to a model introduced by 
Dragulescu and Yakovenko [gj, described by Eqs. 0 with 

A Xj k = ex k - (1 - e)Xj , (9) 

where e is a random number in (0,1]. By comparing 
Eqs. © and ©, one can see that while in the dynam¬ 
ics of the model of Dragulescu and Yakovenko only one 
random number e is present, the model introduced above 
contains two random numbers and reduces to the first 
one by setting ej = 1 — e k - However, in the context of 
wealth exchanges such a strong correlation between ej 
and is difficult to understand. In fact, the dynamics 
of the Dragulescu and Yakovenko model can be better 
understood by rewriting the update relations as 

x'j = e(xj + x k ), 
x k = (1 -e)(xj+x k ), 

representing a random reshuffling in a single time of the 
total initial amount Xj + Xk between the two interacting 
units. 

The presence of two independent random numbers ej 
and Cfc versus the single random number e may seem a 
technical detail, but it implies basically different inter¬ 
pretations. Importantly, the model of Dragulescu and 
Yakovenko leads to the exponential equilibrium wealth 
distribution describing a society where most of the peo¬ 
ple are really poor and the distribution mode is x = 0. 
Instead, the model proposed leads to the T-distribution 
f a (x) with shape parameter a = 2, which corresponds 
to a society where most of the people have the wealth 
around the average value (x) = 1 (assuming that initially 
each unit has a wealth x = 1), with a mode x = 1/2 and 
there is nobody with x = 0. 

The model of Dragulescu and Yakovenko was modified 
by Chakraborti and Chakrabarti [0] assuming that it 
is not the total initial amount Xj + Xk that is reshuffled 
randomly between the two interacting units but only a 
part (1 — A )(xj +Xk), while a fraction A is put aside. The 
corresponding exchange rule reads: 


Xj = A Xj + e(l - A )(xj + Xk ), 
x k = A x k + (1 - e)(l - A )(xj + x k ) • 


( 11 ) 


The latter equations are equivalent to Eqs. 0 with 


Axjk = (1 - A)[(l - e)xj + ex k ] ■ (12) 


In this model the equilibrium distribution of wealth is 
wel l descr ibed by a T-distribution with an a > 1, given 

by " 


1 + 2A 

a = 


1 - A 


(13) 
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if (x) = 1. This implies that f(x = 0) = 0 and a mode 
x > 0, if A > 0. 

Thus, the exponential shape of the equilibrium wealth 
distribution of the Dragulescu and Yakovenko model 
arises from the possibility that during a given interac¬ 
tion a unit can in principle lose all wealth that will go 
to some other unit. Considering Eqs. m, this will hap¬ 
pen to unit j whenever the extracted value of the ran¬ 
dom number e is close enough to zero. The model of 
Chakraborti and Chakrabarti [l3j is in this respect illu¬ 
minating, since it shows that due to the introduction of 
saving (through a saving propensity A > 0) such a sit¬ 
uation never occurs and there are no units for x —> 0, 
leading not to the exponential but to a T-distribution 
(|8|) with shape parameter m Therefore, units car¬ 
rying out immediate exchanges can be formally seen as 
equivalent to units with a saving propensity. More pre¬ 
cisely the immediate exchange model proposed leads to 
the same equilibrium wealth distribution as the model of 
Chakraborti and Chakrabarti with A = 1/4. The fun¬ 
damental difference is that in the model proposed here, 
the situation with Xi = 0 is excluded naturally without 
further assumptions solely by the fact that even if during 
an exchange a unit j gives away everything, he always re¬ 
ceives something in exchange from the other unit k and 
thus, after the transaction, one always finds an Xi > 0 
for each i. 

While the model of Dragulescu and Yakovenko assumes 
that the total amount of wealth, Xj+Xk, is reshuffled ran¬ 
domly between the two interacting units, the first model 
introduced by Angle ||[30| assumes a unidirectional flow: 
with probability po a random fraction e of the wealth Xj 
of unit j is transferred to unit k , while with probability 
1 — Po a random fraction e of the wealth Xk of unit k is 
transferred to unit j. Considering the particular case of 
the symmetrical interaction, i.e., po = 1/2, then the dy¬ 
namical evolution of the model is determined by Eqs. © 
with 




0.01 0.1 1 10 
X 

FIG. 1. Semilogarithmic (top) and logarithmic (bottom) plot 
of the equilibrium distribution p(x) of wealth x for different 
values of p,, i.e., for different fractions of unidirectional inter¬ 
actions. The value g = 0 corresponds to immediate exchanges 
only and g = 1 to the situation when all the interactions are 
unidirectional. Dots represent results from numerical simu¬ 
lations, while lines are a fitting using the T-distribution l|8)l 
with q given by Eq. I© for different values of p. 

with shape parameter [l|] 


Ax jk = e[r)Xj - (1 - rf)x k \ ■ (14) 

Here the stochastic variable g can assume the values 0 
or 1 with probability 1/2; e maintains the same meaning 
as above. In the case of the unidirectional wealth flow 
model of Angle, the growth of the fraction of poor units 
with wealth close to zero is even more dramatic than in 
the reshuffling dynamics considered above: in fact, the 
resulting equilibrium wealth distribution of this model 
is found to be a T-distribution with shape parameter 
a = 1/2 & that diverges for x —>■ 0, signaling a wealth 
accumulation in the hands of a very few units. 

When a saving propensity is introduced in this model, 
assuming that no wealth fraction greater than (1 — A) can 
be transferred from one unit to the other, then Eq. m 
for the exchanged wealth becomes 

Axjk = e(l - X)[vxj - (1 - rf)x k ] ■ (15) 

In this case, the equilibrium wealth distribution for a 
generic value of A is still described by a T-distribution 


+ 2A 
2(“A)’ 


which is just half of the value given by Eq. m for the 
Chakraborti and Chakrabarti model. In this case, the di¬ 
vergence at x —> 0 will persist as long as a < 1, meaning 
saving propensities A < 1 /4. At A = 1/4 the exponential 
distribution is obtained, while for A > 1/4 the distribu¬ 
tion recovers the bell shape with mode x > 0. The value 
a = 2 of the immediate exchange model introduced above 
is recovered for A = 1/2. 

One could expect that (at least in the case of symmet¬ 
rical interactions) a long series of unidirectional wealth 
flows provides results equivalent to those produced by an 
analogous series of immediate (bidirectional) exchanges. 
However, the basic difference in the region at x —> 0 be¬ 
tween the shapes of the equilibrium wealth distributions 
of the model of Angle and of the immediate-exchange 
model proposed here reveal that this is not the case. A 
better understanding of the latter issue can be obtained 
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by investigating a system where at each time iteration 
the pair of units extracted for a trade will carry out with 
probability /z the unidirectional transaction defined by 
Eq. (Ildl) and with probability 1 — /z the immediate ex¬ 
change defined by Eq. ©• 

If n = 0, i.e., all interactions are immediate exchanges, 
the equilibrium wealth distribution is a E-distribution 
with a = 2. For any /z > 0 the equilibrium distributions 
are still well described for sufficiently large x by a E- 
distribution, see Fig. [Q The dependence of a on /z is 
found to be given by 

cr(/z) = 2 exp[—ln(4)/z] = 2 1 ~ 2 ^. (17) 


be exchanged that can vary from time to time, or other 
external random perturbations affecting their decisions. 

Here, the probabilistic acceptance criterion is assumed 
to depend only on the information currently available to 
the two trading units at the moment of the trade, i.e., 
we will use forms of acceptance criterion that depend on 
the wealths owned by the two interacting units before the 
(possible) trade and on the amount of wealth exchanged. 
Without loss of generality, one can focus on the generic 
unit j that will be assumed to accept to carry out the 
trade with an acceptance probability qj. As a simple 
example of acceptance probability we consider in the fol¬ 
lowing the linear piece-wise function given by 


For /z = 1/2, i.e., when half of the interactions are uni¬ 
directional and half are immediate exchanges, a = 1 
and the equilibrium wealth distribution is exponential, 
f(x ) = exp(— x). For x < 0.05 the distributions devi¬ 
ate from the E-distributions, see Fig. [T] for x — 0; the 
deviation is the larger the larger is the fraction of unidi¬ 
rectional interactions, /i. For /z = 1 the model described 
by Eqs. © and (ED is recovered. 

Here, units making immediate exchanges a fraction 
(1 — /z) of times can be interpreted as units with an effec¬ 
tive saving propensity A(/z) < A, where A is the value of 
the saving propensity of the Angle model, corresponding 
to /z = 1. The explicit relation between the effective sav¬ 
ing propensity A(/z) and the parameter /i can be found 
combining Eqs. (llGl) and ED, 


A(/z) 


2 2 ( 1 -/d _ i 
2 2 (i-/d + 2 ' 


(18) 


From here one can see that to /z = 1 corresponds A = 0, 
i.e., the units can lose everything during an exchange. In¬ 
stead, /z = 1/2 leads to the same equilibrium distribution 
as the model of Angle with saving propensity A = 1/4, 
and /z = 0 to the equilibrium distribution with saving 
propensity A = 1/2. Distributions corresponding to val¬ 
ues of A > 1/2 cannot be obtained, since they would 
correspond to negative values of /z. 


III. THE INFLUENCE OF A TRADING 
CRITERION 

A. Formulation of the acceptance criterion 

Let us now introduce a trading criterion on the basis 
of what each unit will decide whether to make the trade 
or not (in other kinetic wealth exchange models the two 
randomly chosen units always make the trade). The in¬ 
troduction of probabilistic factors influencing trades is 
a possible way to go beyond the assumption of perfect 
knowledge of the units assumed in neo-classical economic 
models. A probability law suitably describes the natural 
lack of perfect knowledge of the trading units concerning 
the product and its actual value measured as wealth, as 
well as the effect of personal feelings about the goods to 


{ 0, if Ax jk < —rj(x ), 

1 + A x jk /(r)(x)) , if - r/(x) < Ax jk < 0, 
1, if A Xj k > 0 . 

(19) 

Here 77 is a parameter, (x) = ]TL = 0 X i/N is the average 
wealth of the system, and A Xj k is given by Eq. ©. This 
function is depicted as a continuous line in Fig. [2] 



FIG. 2. Examples of probability qj(Axjk) that unit j accepts 
to make a wealth exchange with unit k leading to a net wealth 
variation A Xjk of unit j. Here “linear” refers to the piecewise 
linear probability in Eq. ED for 7] = 1 / 2 , while the other 
“exponential” lines represent the probability given by Eq. ED 
for rj = 1/2 and 5xo = —0.2; 0; 0.2. See text for details. 

In this example, the criterion for accepting or not the 
exchange depends on the quantity r/(x) and on the ex¬ 
changed wealth Axjk- If A Xj k is negative (the unit j 
is going to lose something) with | Aoryfc | comparable with 
the scale ?/(lc), then it is highly probable that there will 
be no trade as the unit j who would lose is most likely 
not interested of the exchange. The exchange will cer¬ 
tainly not take place if |Acnyfc | > r](x). If Axj k is still 
negative but with \Axj k | not too large compared to the 
scale rj(x) then it is possible that unit j will accept the 
exchange anyway, even though losing something. The 
latter situation could correspond, e.g., to a situation in 
which he is not capable to estimate the loss, or he is 
simply interested in the other object for whatever rea¬ 
son and therefore is ready to accept a limited loss. For 
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A Xjk > 0 unit j will always accept the exchange. Both 
trading units j and k will independently adopt the same 
criterion, the criterion of unit k being defined by an anal¬ 
ogous function q k (Axkj) = qj(—Axjk) obtained from the 
function given by Eq. (fl9l) by exchanging k and j and, 
correspondingly, Axkj with —Axjk- The decision of each 
agent k or j is taken probabilistically with the help of 
an additional random number that is extracted indepen¬ 
dently and compared with qj(Axjk) or q k {Axkj) for the 
decisions of unit j and k, respectively. Since the quan¬ 
tity A Xjk = CkXk — ejXj is given as the difference between 
the values e k Xk and ejXj of the wealth goods exchanged, 
it represents an actual measure of the profit or loss of 
a unit in an exchange. Thus, the type of interaction 
between units formulated above is expected to simulate 
an actual barter or trade realistically, since whether the 
units will carry out the exchange or not depends on their 
corresponding profit or loss. 

It is possible to construct similar functions with spe¬ 
cific properties and depending on additional parameters. 
We employ in the following some other forms of the ac¬ 
ceptance probability function, in particular the following 
exponential shape, 



0 12 3 4 


x 

FIG. 3. Equilibrium wealth distributions f(x) of the 
immediate-exchange model. All the distributions obtained 
using different or no acceptance criteria (dots) collapse on 
the same E-distribution f a (x) with a = 2 (continuous 
curve). The criteria used employ (a) the “linear” proba¬ 
bility in Eq. (1191) with q = 0.1; 0.5; 1; 5; 10; and (b) the 
“exponential” probability in Eq. HU) with q = 0.5 and 
Sx o = —0.4; —0.2; 0; 0.2; 0.4. Inset: same plot in semiloga- 
rithmic scale. 


qj (A Xjk) 


exp [(Ax - 8x 0 )/{q{x ))\, if Ax jk < Sx 0 , 
1, if A Xjk > Sx o . 


( 20 ) 

This probability function has an exponential shape for 
A Xjk < Sx o and is equal to 1 otherwise. The shapes 
of the curve corresponding to an acceptance parameter 
q = 0.5 and shift parameters Sxq = 0, ±0.2 are shown 
in Fig. [2] The parameter Sx o can be used to make the 
acceptance criterion stricter or looser - notice that an 
analogous parameter can in principle be introduced in 
any other function qj(Axjk)- The acceptance criterion 
become looser for Sxo < 0, when a unit can accept to 
carry out an exchange even if accompanied by a limited 
loss of the order of |feo|. Instead, the criterion becomes 
stricter for Sx o > 0, since in this case trades are accepted 
with certainty only if the gain is larger than the threshold 
value fap. 


B. Acceptance criterion in the immediate exchange 

model 

We start by applying the acceptance criterion m to 
the kinetic model of immediate exchange introduced in 
Sect. HD The equilibrium distributions of wealth obtained 
after adding the acceptance criterion turns out to be a 
E-distribution ([5]), namely the distribution / 2 (x) with 
the same shape parameter a = 2 independently of the 
value of q. Thus, the shape of the equilibrium wealth 
distribution remains unchanged with respect to the case 
in which no acceptance criterion is used. 

This is true also when different acceptance criteria are 
employed. The equilibrium wealth distributions obtained 
employing the piece-wise linear function (1191) and those 


corresponding to the exponential acceptance probability 
function (l20l) for different values of the shift parameter 
Sx o are compared with each other and with the analytical 
shape of the distribution / 2 (x) in Fig. [3] 

The value of q in Eqs. m and (USD as well that of 
Sx o in Eq. (l20l) change the relaxation process to equi¬ 
librium [32[ (e.g. a smaller q corresponds to a larger 
relaxation time) but surprisingly do not have any impor¬ 
tance for the equilibrium wealth distribution, that is ex¬ 
actly the same one obtained when assuming that the two 
randomly chosen units make the trade with probability 

q= I- 


C. Asymmetrical acceptance criteria 

The considerations above lead to the conclusion that 
introducing a decision making process that can be inter¬ 
preted as the trial to introduce intelligence in the units 
behavior, is irrelevant for the final state of the system. 
However, this is true only as long as the acceptance crite¬ 
rion is formulated symmetrically with respect to the two 
trading units. Different asymmetrical criteria lead to dif¬ 
ferent forms of the equilibrium wealth distribution and 
to shapes of the equilibrium wealth distributions different 
from E-distributions. Such asymmetrical criteria can be 
of interest in the study of preferential attachment effects, 
e.g., when considering the influence of the richness of the 
two trading units on the outcome of the trade. That is 
the case of the so-called “rich get richer effect” already 
considered in the model of Angle Q, characterized by 
uni-directional wealth flows taking place with a higher 
probability from the poorer to the richer unit. Notice 
that unidirectional models represent actual realizations 
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of the so-called “Matthew effect” (33|, since during an 
encounter not only the richer unit becomes richer with 
a higher probability but correspondingly, due to wealth 
conservation, the poorer unit becomes even poorer. 

As a test, we have studied an asymmetrical version of 
the immediate exchange model discussed above, in which 
a criterion favoring richer units was introduced by always 
allowing trades with a gain for the richer trading unit 
and preventing trades accompanied by a net gain of the 
poorer unit a fraction 9 of times, 0 < 9 < 1. 

In Fig. U the equilibrium wealth distribution for some 
values of 9 ranging between 9 = 0 (corresponding to sym¬ 
metrical exchanges) and 9 = 0.9 (representing a strong 
asymmetry in the exchange criterion) are shown. One 
can observe how the number of units in the small-a; re¬ 
gion increases with 9 and that, correspondingly, the dis¬ 
tribution mode shifts leftwards. It is to be noticed that 
for high enough values of 9 the distribution diverges for 
x —> 0 — Fig. 2] shows such a case for 9 = 0.9. Most 
importantly, the shapes of the equilibrium distributions 
obtained for 9 > 0 are not well fitted anymore by a f- 
distribution (comparison not shown). 



0 12 3 4 5 


x 

FIG. 4. Equilibrium wealth distribution obtained using an 
asymmetrical criterion for different fractions 9 of trades that 
would be advantageous for the richer unit; see text for details. 
Equilibrium distributions with 9 > 0 are not well fitted by a 
F-distribution (comparison not shown). 


D. Acceptance criteria in other models: relative 
acceptance criteria 

In order to check for other instances of the invari¬ 
ance properties found above for the immediate-exchange 
model, we have numerically investigated the effects of 
an acceptance criterion in other kinetic exchange mod¬ 
els, namely in the Dragulescu-Yakovenko Q and in 
the Chakraborti-Chakrabarti model with saving propen¬ 
sity 0 . 

We start from the Dragulescu-Yakovenko model, de¬ 
fined by Eqs. (El or by Eqs. © with ©, whose equilib¬ 
rium wealth distribution is a simple exponential shape, 


f(x) = exp(-x) (if (x) = 1 ) In fact, we find that 
the exponential form remains unchanged, independently 
of the form of the symmetrical acceptance criterion used, 
such as those defined by Eqs. m or CUl) . The indepen¬ 
dence from the acceptance criteria widens the validity 
range and justifies the use of the original minimal version 
of the model where no acceptance criterion is present, in¬ 
stead of more elaborate versions where units follow more 
“intelligent” principles. 

However, the independence of the equilibrium distri¬ 
bution of a symmetrical acceptance criterion does not 
hold for the Chakraborti-Chakrabarti model with saving 
propensity fl3 ] . In this case one obtains a different equi¬ 
librium wealth distribution, that cannot be fitted any¬ 
more by a F-distribution. 

Considering the homogeneity of the update rules of 
the model, we have checked the effect of a probabilis¬ 
tic constraint on the relative variable A Xjk/xj of unit j 
rather than on the absolute wealth exchange A Xjk, find¬ 
ing that a good fitting of the equilibrium wealth distribu¬ 
tions through a F-distribution is recovered in this case. 
In other words, if in the Chakraborti-Chakrabarti model 
with a given saving propensity A an acceptance criterion 
is introduced, defined in terms of an (arbitrary) accep¬ 
tance probability function Qj for each unit j depending 
on the relative amount of wealth exchanged A Xjk/xj, 
then the equilibrium wealth distribution is still a F- 
function, but with a shape parameter different from the 
value given by Eq. m■ For clarity we discuss only the 
acceptance criterion defined by the acceptance probabil¬ 
ity Qj obtained by modifying the piece-wise linear accep¬ 
tance probability qj in Eq. (fliil) . 


Qj (A Xjk/xj ) — < 


0 , 

1 + A Xj k /{rj Xj), 


if A Xjk/xj < - 77 , 
if — rj < A Xjk/xj < 0, 


1, if A Xjk > 0 . 

( 21 ) 

The rescaled variable A Xjk/xj has the new range (—1,1). 
Results are summarized in Fig. [SJtop, showing the new 
shape parameter a q of the F-distribution f aq (x) fit¬ 
ting the equilibrium wealth distribution versus the corre¬ 
sponding value of the saving propensity A of the model, 
for different values of the acceptance parameter 77 . In the 
limit of very large values of 77 , in which the acceptance 
criterion become very loose and becomes equivalent to 
no criterion present, results reduce to those of the orig¬ 
inal Chakraborti-Chakrabarti model, represented by the 
continuous line representing a(X) of Eq. (fl3l) ; in all other 
cases one has a q > a(A). 


The results can also be expressed by saying that the 
relative acceptance criterion with acceptance parameter 
77 , defined by Eqs. m, turns the equilibrium wealth 
distribution of the Chakraborti-Chakrabarti model with 
saving propensity A into that of the same model with a 
different saving propensity X q . The value of the effective 
saving propensity X q corresponding to the observed a q 
can be obtained inverting the relation (1131) and is given 









by 


100 


A, = (a,-l)/(a, + 2). (22) 

For convenience the values of X q are depicted in Fig. [5} 
bottom. It is clear that for every value of the acceptance 
parameter 77 one has X q > A, i.e., the use of an accep¬ 
tance criterion is equivalent to a larger saving propen¬ 
sity. For large values of 77 , when the acceptance crite¬ 
rion becomes very loose, results reduce to the original 
Chakraborti-Chakrabarti model without acceptance cri¬ 
terion and X q —> A. 

A surprising point implied by these results is that 
the saving propensity of the Chakraborti-Chakrabarti 
model can be generated by introducing a relative accep¬ 
tance criterion. In fact, the considerations above also 
apply to the Dragulescu-Yakovenko model, that is ob¬ 
tained from the Chakraborti-Chakrabarti model in the 
limit A —> 0. Thus, one can start from the Dragulescu- 
Yakovenko model without saving propensity and obtain 
the same equilibrium distribution of the Chakraborti- 
Chakrabarti model with arbitrary saving propensity A, 
corresponding to the dots at A = 0 in Fig.jSj through the 
use of an acceptance criterion with a suitable value of 
77 . This justifies at the microscopic level of the decision 
processes, carried out by the units before an exchange, 
the use of a saving propensity A > 0 inserted “by hand” 
in the Chakraborti-Chakrabarti model as a convenient 
numerical or effective statistical procedure. 

IV. HETEROGENEOUS 
IMMEDIATE-EXCHANGE MODEL 

In this section we investigate the generalization of the 
homogeneous model studied above to the heterogeneous 
case, when each unit behaves according to a different 
trading criterion. As in the homogeneous case, the pres¬ 
ence of a heterogeneity in the parameter values can be 
related to specific features of the units, such as the effi¬ 
ciency of the trading strategy employed by the units or 
their education level 

Here, we limit ourselves to consider an absolute accep¬ 
tance criterion. In fact, while the equilibrium distribu¬ 
tion of a homogeneous system is invariant under change 
of the form of the (symmetrical) criterion employed, the 
presence of a heterogeneity level in the system parame¬ 
ters is sufficient to break this invariance. The same piece- 
wise linear probability function considered above is used, 
but now each unit j (j = 1 ,N) is characterized by 
a different parameter rjj. The corresponding acceptance 
probability function is given by 

{ 0, ii Ax jk <-r]j(x), 

1 + Axjk/(r}j(x)) , if - r}j{x) < Axj k < 0, 

1, if Axjk > 0. 

(23) 

The equilibrium distribution does not have a universal 
shape anymore and depends on the specific acceptance 




FIG. 5. Fitting of the equilibrium wealth distribution of the 
Chakraborti-Chakrabarti model with acceptance criterion de¬ 
fined by Eq. m through a E-distribution. Top: shape pa¬ 
rameter a q versus saving parameter A for different acceptance 
parameters 77 . The black solid line is the value of a(A) given 
by Eq. (O, corresponding to the Chakraborti-Chakrabarti 
model, recovered in the r/ —»• 00 limit, i.e., when there is no 
acceptance criterion. Bottom: effective saving parameter X q 
corresponding to the values of a q in the top panel. In general 
X q > X. In the limit of large 77 the presence of an acceptance 
criterion becomes negligible and X q — > X. 


parameters { 77 ^} used or on their distribution (j>{rf) that 
can be defined in the limit of a large number of units. 

It to be noticed that heterogeneity does not concern 
the functional form of the amount of exchanged wealth 
Axjk but enters the model through the parameters of the 
acceptance criterion used by each unit. 

We have investigated various forms of distributions 
4>{r]). The equilibrium wealth distributions obtained from 
uniform threshold distributions with parameters rji ran¬ 
domly extracted in an interval (77 m j n , rj max ) in general are 
not realistic. An example of such equilibrium wealth dis¬ 
tributions obtained for a fixed 77 m j n = 0.1 and different 
values of rj max are plotted in Fig. [Gj One can notice that 
when thresholds are very limited in the range (all close to 
77 = 0.1 in the example of Fig. [ 6 ]) the model is almost ho¬ 
mogeneous and one recovers a E-distribution. The pres¬ 
ence of a significant fraction of agents i with higher values 
of the parameter r]i produces a depletion of the distribu- 
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tion in the intermediate wealth zone and an enhancement 
of the concentration of units both at small and large val¬ 
ues of wealth. In the large-wealth region one can notice 
the formation of a new maximum that can be interpreted 
as the formation of a rich class made of the units with the 
smallest values of r/i’s. To the best of our knowledge this 
type of maximum is not observed in real data of wealth 
distributions. 

Results reproducing a realistic shape of wealth distri¬ 
butions, are obtained by following a simple recipe con¬ 
sisting in constructing a system in which the majority of 
units has a standard or low performance during trading 
interactions and a small fraction of units perform much 
better. This prescription was already tested in Ref. [HJ 
in the framework of the Chakraborti-Chakrabarti model 
with saving propensities [l3| by assigning a zero saving 
propensity to 99% of the population and higher saving 
propensities, uniformly distributed in A S (0,1), to the 
remaining 1% of the population. This led to realistic 
shapes of wealth distribution, see Ref. 0 for details. 
The recipe is here adapted to the immediate-exchange 
model under consideration by assigning a (fixed) larger 
parameter value r/i = fj = 2 to a majority (95%) of the 
population here a larger value of rj corresponds to 
a looser criterion used in accepting an exchange, while 
the remaining 5% of the population is assigned a set of 
smaller parameters r/i uniformly extracted in the inter¬ 
val r/ G (0.5, 0.7). The resulting equilibrium wealth dis¬ 
tribution, shown in Fig. [71 indeed presents the stylized 
features of wealth distributions, such as a mode x m > 0 
and a Pareto power-law at larger values of wealth with a 
Pareto index p « 2. 



X 

FIG. 6. Equilibrium distribution f(x) of wealth a; of a system 
of heterogeneous units with diversified acceptance parame¬ 
ters {r/i} uniformly distributed in the interval (r/min, r/max), 
for r/min = 0.1 and different values of r/ max • 


V. CONCLUSIONS 

In the present paper we have proposed a model of 
wealth exchange aimed at overcoming some critics raised 



FIG. 7. Equilibrium distribution f(x) of wealth a: of a system 
of heterogeneous units with diversified acceptance parameters 
{r/i}'- 95% of the units has r/ — 2, while the remaining 5% has 
{r/i} uniformly distributed in r/ € (0.5, 0.7). 


about the kinetic exchange models, laying down firmer 
foundations for the models themselves, and providing a 
more realistic description of a trade between two eco¬ 
nomic units. Doing this, two main features of the orig¬ 
inal models have been modified. First, an immediate- 
exchange wealth dynamics has been introduced, describ¬ 
ing two actual and distinct wealth flows taking place dur¬ 
ing the same interaction between two economic units, just 
as in an actual trade. Models, where the total wealth 
is randomly reshuffled between the two units or which 
follow a uni-directional dynamics, better represent ex¬ 
changes of valuables without the demand of an immedi¬ 
ate reward, as, e.g., in gift economies. Immediate bidirec¬ 
tional exchanges produce more realistic shapes of wealth 
distribution /(a;), such that /(0) =0 and having a mode 
larger than zero. Furthermore, the model is made more 
realistic at a microeconomic level, allowing units to de¬ 
cide whether to trade or not, according to some criterion 
based on the profit gained from the trade. Criteria for¬ 
mulated symmetrically with respect to the two trading 
units — but arbitrary for the rest — preserve some sym¬ 
metries of the model. In the immediate-exchange models 
as well as in the Dragulescu-Yakovenko model, a sym¬ 
metrical criterion does not change the equilibrium wealth 
distribution at all. 

In the case of the Chakraborti-Chakrabarti model, 
when a symmetrical acceptance condition on the rela¬ 
tive wealth is used, one still finds that the equilibrium 
distribution is a T-distribution, but with a larger value 
of the shape parameter a, corresponding to a higher sav¬ 
ing propensity A. This implies interesting links between 
different kinetic exchange models and provides at the 
same time a microscopic justification. For example, the 
equilibrium distribution of the Chakraborti-Chakrabarti 
model can also be obtained from the reshuffling dynamics 
of the Dragulescu-Yakovenko model with the addition of 
an acceptance criterion with suitable parameters. 

The situation changes either when the criterion de- 
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pends asymmetrically on the wealths of the two units 
(also in homogeneous systems) or in the case of hetero¬ 
geneous systems, where there is anyway a situation of 
asymmetry due to the fact that each unit uses a dif¬ 
ferent criterion. In these cases the equilibrium wealth 
distribution depends on the specific set of criteria used 
by units. The heterogeneous version of the bidirectional 
model presents a wide spectrum of possible shapes for 
the equilibrium wealth distribution. We have found that 
realistic shapes, presenting a Pareto power-law, are ob¬ 
tained for moderate levels of heterogeneity. 

In order to provide quantitative explanations of the 
many regularities of and links between the various models 
discussed above it would be important that future work 
advances also toward the analytical solution of kinetic ex¬ 
change models. Furthermore, while the ubiquitous pres¬ 


ence of the F-distribution can be related to the validity 
of the Boltzmann theorem, due to the conservation of 
wealth during trades, it would be relevant to gain some 
information also on the analytical shapes of the equilib¬ 
rium wealth distributions when the Boltzmann theorem 
does not apply, e.g . in the interesting Maxwell-demon¬ 
type dynamics [29( corresponding to cases of asymmet¬ 
rical acceptance criteria describing e.g. the “rich gets 
richer” effect. 
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